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[i]  Under-dense  HF  heating  experiments  were  conducted  near  local  solar  noon  as  well 
as  in  the  nighttime  with  the  HF  heater  transmitting  at  9. 1  MHz  directed  along  the 
geomagnetic  zenith  and  run  at  2  min  on  and  2  min  off.  The  effective  isotropic  radiated 
power  of  the  HF  transmitter  exceeded  3  GW.  The  Digisonde  operated  in  a  fast  mode 
was  used  to  monitor  the  temporal  evolution  of  the  ionospheric  electron  density 
distributions  in  the  bottomside  of  the  ionosphere  (in  the  ranges  from  90  to  190  km  in 
the  noontime  and  from  230  to  350  km  in  the  nighttime).  The  electron  temperature 
distributions  were  then  evaluated.  The  results  show'  that  the  electron  density  distributions 
arc  modified  continuously  over  the  experimental  periods.  In  the  noontime,  the  electron 
density  deercases/inereases  in  time  in  the  region  below/above  a  height  at  about 
140  km,  manifesting  the  ehange  of  the  balance  between  the  photoionization  and  the 
electron-ion  recombination  and  the  electron-oxygen  dissociative  attachment  losses  by 
the  heating.  In  the  nighttime,  the  ionosphere  was  lifted  by  30  to  50  km  through 
continuously  upward  expansion,  resulting  in  the  drop  of  the  electron  density  in  the 
bottomside  of  the  ionosphere  in  time.  A  comparison  w  ith  the  ionogram,  height,  and 
electron  density  distribution  of  unheated  ionosphere  with  similar  background  conditions 
further  elaborates  the  observation  of  thermal  expansion. 
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1.  Introduction 

[2]  The  Ionosphcnc  Research  Instrument  (IRI)  of  the  High 
Frequency  Active  Auroral  Research  Program  (HAARP)  is 
a  HF  transmitter,  which  delivers  0.36  to  3.6  GW  effective 
isotropic  radiated  powers  (EIRP)  for  the  radiation  frequen¬ 
cies  from  2.8  to  10  MHz.  The  power  density  of  the  HF  wave 
in  the  ionosphere  is  up  to  about  30  mW/m2  at  100  km  alti¬ 
tude.  It  is  anticipated  that  the  HF  heating  of  the  ionosphere 
via  both  collision  and  eollisionless  heating  processes  will 
make  a  significant  impact  on  the  plasma  dynamics  of  the 
ionosphere  [ Fejer ,  1979;  Gurevich ,  2007;  Kuo  el  al. ,  2009; 
Pedersen  el  al ,  2009;  Kuo  el  al .,  2010;  Kuo  and  Snyder , 
2010],  which  provides  useful  information  for  studying  the 
properties  and  behavior  of  the  ionosphere. 

[3]  In  the  past,  most  of  the  HF  heating  experiments  were 
focused  on  overdense  situation  by  applying  O  mode  HF 
heaters  with  frequencies  less  than  the  foF2.  In  this  situation, 
parametric  instabilities  can  be  excited  locally  near  or  below 
the  HF  reflection  height  to  contribute  to  anomalous  heating, 
where  the  frequency  and  wave  vector  matching  conditions 
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for  parametric  coupling  can  be  satisfied  [Fejer,  1979;  Stenjlo, 
1985;  Kuo,  2001].  Parametric  instabilities  have  been  respon¬ 
sible  for  the  observation  of  new  physical  phenomena  in  the 
experiments  [Gordon  and  Carlson,  1974;  Stuhbe  el  al.,  1984, 
1992].  Moreover,  the  anomalous  heating  expanding  along 
the  geomagnetic  field  may  cause  ionospheric  modification 
over  a  large  altitude  region.  In  particular,  by  combining  solar 
illumination  in  the  solar  noon  time  period  with  the  anomalous 
heating  of  the  HF  heater,  significant  modification  of  the 
ionospheric  plasma  density  distribution  was  observed  [Kuo  and 
Snyder,  2010],  which  could  impact  communication  systems. 

[4]  In  the  case  of  under-dense  heating,  the  EIRP  of  the 
HAARP  heater  can  be  increased  significantly  by  increasing 
the  heater  frequency  With  higher  heater  frequency,  the  loss 
of  heating  power  in  the  D  region,  in  the  daytime  experiment, 
is  also  reduced.  In  the  daytime,  it  is  anticipated  that  the 
electron  heating  will  reduce  the  clcctron-ion  recombination 
coefficient  to  increase  the  plasma  density,  provided  that 
transport  is  insignificant  However,  in  the  E  region  the  rate 
of  attachment  of  electrons  to  oxygen  can  also  be  increased, 
resulting  in  decreasing  plasma  density  In  the  nighttime,  the 
ionosphere  may  be  modified  via  thermal  expansion.  There¬ 
fore,  the  modification  of  the  ionosphere  via  overdense 
and  under-dense  HF  heating  as  w'ell  as  in  the  daytime  and 
nighttime  can  be  quite  different 

[5]  This  work  reports  experiments  conducted  with  the 
9.1  MHz  HF  heater  that  explores  the  under-dense  heating 
effects  on  the  ionospheric  electron  density  distributions  in 
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Figure  1.  Comparison  of  common  (to  on/off)  electron  density  distributions  with  consecutively  increas¬ 
ing  heating  times.  The  one  corresponding  to  the  first  on/off  period  from  2032  to  2036  UTC  is  used 
as  a  reference. 


the  solar  noon  time  period  and  in  the  nighttime  for  a  com¬ 
parison.  The  inelastic  collision  processes  are  revealed  via 
the  change  of  the  electron  density  distribution  in  the  noon¬ 
time  heating  experiment. 

2.  Experiment 

[6]  On  12  April  2010  from  2032  to  2104  UTC  (1232  to 
1304  local  time)  and  on  13  April  from  0812  to  0844  UTC 
(0012  to  0044  local  time),  using  the  HAARP  transmitter 
facility  at  Gakona,  AK,  at  full  power  (3.6  MW),  experiments 
were  conducted  with  the  HF  heater  transmitting  at  9. 1  MHz 
directed  along  the  geomagnetic  zenith  and  run  at  2  min  on 
and  2  min  off.  In  the  on  periods,  the  polarization  of  the 
heating  wave  was  switched  alternately  with  O  mode  and  X 
mode.  The  12*15  array  of  the  transmitter  was  split  into  two 


6x15  subarrays;  one  was  run  at  ew  full  power  and  the  other 
run  at  100%  power  modulation  (PM)  by  rectangular  waves 
of  2.02  kHz,  5  kHz,  8  kHz,  and  13  kHz.  Each  modulation 
frequency  was  applied  consecutively  in  two  2  min  on  per¬ 
iods  at  the  sequence  of  O/X  mode  heater  and  of  100%  power 
modulation  to  one  of  the  subarray.  The  purpose  of  modu¬ 
lation  is  to  generate  VLF  wave  by  wave  beating  through 
the  nonlinearity  of  plasma.  While  the  unexpected  high-level 
background  noise  at  13  kHz  contaminated  VLF  data,  new 
phenomena  directly  related  to  HF  heating  of  the  background 
plasma  were  observed.  These  new  phenomena  do  not  depend 
on  the  polarization  of  the  HF  heater  because  the  heater 
frequency  is  much  larger  than  the  peak  plasma  frequency 
(foF2).  The  experiment  was  monitored  by  the  HAARP 
Digisonde.  The  results  developed  from  the  ionogram  data 
files  are  presented. 
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Figure  2.  Time  evolution  of  the  electron  density  at  three 
heights:  1 10,  140,  and  170  km  in  the  noontime  experiment 

[^]  The  Digisondc  acquired  ionograms  every  minute  dur¬ 
ing  the  heater  on  and  off  periods.  The  1.0  to  5.5  MHz 
ionograms,  with  75  kHz  steps,  started  on  the  minute.  The 
recording  of  each  lonogram  took  less  than  19  s  Thus,  the 
ionogram  acquisition  time  from  approximately  2  to  4.6  MHz, 
the  frequency  band  of  interest,  was  less  than  1 1  s. 

3.  Experimental  Results 
3.1.  Noontime  Experiment 

[s]  Quantitative  frequency  distribution  and  time  variations 
of  the  virtual  height  and  height  spread  of  the  sounding 
echoes  can  be  examined  from  the  ionogram  data  files.  It  was 
found  that  during  the  experiment,  the  virtual  height  traces  of 
the  echoes  in  the  ionograms  were  modified  continuously. 
HF  heater  changed  the  virtual  height  and  spread  of  the 
sounding  echoes.  In  the  analysis,  we  choose  data  only  from 
the  vertical  O  mode  sounding  echoes.  Wc  divide  the  O 
mode  sounding  echoes  in  each  pair  of  heater  on/off  iono¬ 
grams  into  three  groups.  These  three  groups  arc  (1)  unique 
off  ionogram  echoes,  which  appear  only  in  the  heater  off 
ionograms;  (2)  unique  on  ionogram  echoes,  which  appear 
only  in  the  heater-on  ionograms;  and  (3)  common  ionogram 
echoes,  which  appear  in  both  the  hcater-on/off  ionograms 
being  compared.  The  unique  on/off  echoes  are  the  additional 
spread  echoes  in  virtual  height  to  those  of  the  common 
echoes  in  the  heater-on/off  periods.  Their  difference  man¬ 
ifests  the  heater  induced  change  of  density  irregularities 
which  are  present  in  the  background  plasma.  The  virtual 
height  of  the  strongest  received  signal  at  each  frequency 
is  then  taken  to  obtain  the  frequency  distribution  of  the 
virtual  height  of  the  common  echoes  in  each  on-off  iono¬ 
gram  pair. 

[9]  The  virtual  height  distribution  of  the  O  mode  sounding 
echo  does  not  represent  the  actual  height  distribution  of 
the  ionospheric  plasma.  The  University  of  Massachusetts 
developed  Digisondc  provides  a  software  program  SAO 
Explorer  [Galkin  et  a/.,  2008]  (see  http://ulcar.uml.edu/ 
SAO-X/SAO-X,html)  for  rescaling  of  the  virtual  height 


trace  h'(f)  and  calculation  of  the  true  height  profile  h(J).  This 
program  is  used  to  convert  a  virtual  height  distribution  of 
common  ionogram  echoes  to  a  true  height  distribution 
versus  frequency,  which  is  then  inverted  to  the  dependence 
of  electron  density  on  altitude. 

[10]  The  common  electron  density  distributions  in  the 
8  on/off  periods  are  evaluated  from  the  data  files  of  common 
ionogram  echoes.  The  one  corresponding  to  the  first  oatoff 
period  from  2032  to  2036  UTC  was  used  as  a  reference 
for  showing  the  temporal  evolution  of  the  electron  density 
distribution  in  the  subsequent  on/off  periods.  As  shown  in 
Figure  1,  the  common  electron  density  distribution  evolves 
continuously  in  time.  It  suggests  that  the  heater-induced 
modification  of  the  electron  density  distribution  lasts  longer 
than  the  2  min  off  period,  and  the  modification  accumulates 
in  time.  It  also  shows  that  after  a  few  on/off  periods,  the 
density  profile  reaches  a  stable  feature,  the  electron  density 
becomes  lower/higher  than  that  of  the  reference  in  the 
lower-/highcr-altitude  regions 

[11]  The  temporal  evolution  of  the  electron  densities  at 
three  representative  altitudes:  1 10,  140,  and  170  km,  is  pre¬ 
sented  in  Figure  2  to  showr  the  contrast  of  density  changes 
in  time.  As  shown,  the  density  at  110/170  km  decreases/ 
increases  with  time  after  a  transient  period  and  the  density  at 
140  km  varies  slightly  around  a  constant  level.  The  unique 
on/off  echoes  are  the  additional  echo  spreads  from  the 
common  echoes  in  the  on/off  situations  and  their  difference 
manifests  the  heating  effect.  Thus  the  frequency  dependen¬ 
cies  of  the  virtual  heights  of  the  unique  oa^off  echoes  are  also 
converted  to  the  frequency  dependencies  of  the  true  heights 
of  the  unique  on/off  echoes,  which  determine  the  electron 
densities  at  the  true  heights  of  the  unique  oa/off  echoes. 
Combining  the  electron  densities  at  the  true  heights  of  the 
common  echoes  and  of  the  unique  on/off  echoes,  wc  have 
generated  “modeled”  electron  density  distributions  at  heater 
on/off,  which  contain  spatial  density  variations  giving 
“equivalent  effects”  on  the  virtual  height  spreads  of  the 
sounding  echoes.  Normally,  one  can  choose  a  single  virtual 
height  for  each  strongest  received  signal  at  each  frequency  to 
obtain  a  frequency  distribution  of  the  virtual  height  for  the 
unique  on  or  unique  off  echoes,  however,  in  a  few  cases 
observed  in  the  present  experiment,  the  virtual  height  of  the 
strongest  received  signal  has  double  values  with  large  height 
separations.  This  is  demonstrated  in  Figure  3a,  which  sug¬ 
gests  that  the  lonosondc  can  trace  more  than  a  single  line 
when  the  ionospheric  density  distribution  does  not  have  a 
sharp  horizontal  layer  due  to  the  presence  of  field-aligned 
density  irregularities.  A  set  of  true  height  results,  for  the  on/ 
off  at  2041/2043  UT,  based  on  the  virtual  height  distributions 
presented  in  Figure  3a  is  presented  in  Figure  3b  As  shown, 
the  heater  has  modified  the  magnitude  and  the  scale  length  of 
the  effective  electron  density  variations  of  the  modeled 
electron  density  distributions,  which  are  in  quasipcriodic 
forms  resembling  density  irregularities.  Note  that  each  line  in 
Figure  3b  intersects  with  a  horizontal  line  only  once;  thus, 
these  plots  show  only  density  variations  in  altitude  (i.c.,  in 
vertical  direction).  The  scale  length  L  of  the  density  varia¬ 
tion  in  vertical  direction  can  be  used  to  determine  the  scale 
length  Lg  of  the  field-aligned  density  irregularities  by  the 
relation  Lg  =  L  cos0d,  where  0d  is  the  magnetic  dip  angle  and 
equals  76°  at  Gakona,  Alaska.  The  present  approach  models 
the  heating-induced  ionospheric  modification  manifested  by 
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Figure  3.  (a)  Virtual  heights  of  the  strongest  unique  on/off  and  common  echoes  versus  the  sounding 
frequency  and  (b)  comparison  of  the  electron  density  distributions  in  heater  on/off  periods,  including 
the  common  part  as  a  reference. 


the  virtual  height  spread  of  the  sounding  signal  via  an 
equivalent  perturbation  on  the  density  distribution. 

3.2.  Nighttime  Experiment 

[12]  The  bottomsidc  electron  density  distributions  in  the 
8  on/off  periods  from  0812  to  0844  UTC  on  13  April  are 
evaluated  in  the  same  way.  Four  of  the  16  plots  in  the  X 
mode  heater  on  periods  are  presented  in  Figure  4a  to  show' 
the  temporal  evolution  of  the  electron  density'  distribution. 
As  shown  the  distribution  shifts  up  consecutively  with 
decreasing  density  as  the  total  heating  time  increases.  The 
temporal  evolutions  of  electron  densities  at  five  heights: 
230,  255,  280,  305,  and  330  km  are  presented  in  Figure  4b. 
In  the  initial  on/ofFon  period,  the  density  dccrcascs/in- 
ereascs  in  the  lower/highcr  altitude;  then  the  trend  of  density 


decreasing  in  time  after  this  initial  period  is  clearly  shown 
The  transient  period  covering  the  first  oa^ff/on  suggests 
that  it  takes  up  to  6  min  for  the  ionosphere  to  reach  a  steady 
thermal  expansion,  which  upshifts  the  ionosphere  and  con¬ 
tinuously  depletes  the  bottomsidc  plasma  density.  Figure  4b 
also  shows  that  the  densities  drop  quickly  in  the  beginning, 
e.g.,  from  0818  to  0830  UT,  and  those  in  the  lower  altitude 
reach  steady  state  levels  first.  This  trend  is  consistent  with 
that  in  an  upward  thermally  expanding  ionosphere 

4.  Physical  Interpretation 

[13]  The  HF  heater  of  9.1  MHz  cannot  excite  parametric 
instabilities  for  anomalous  heating,  the  wave  heating  was 
limited  to  eollisional  heating.  The  electron  density  distribution 


Figure  4.  (a)  Comparison  of  electron  density  distributions  with  a  consecutive  increase  of  the  total  heat¬ 
ing  time  in  the  nighttime  experiment  and  (b)  time  evolution  of  the  electron  density  at  five  heights:  230, 
255,  280,  305,  and  330  km. 
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Electron  Temperature  (K) 

Figure  5.  Comparison  of  the  unperturbed  electron  temper¬ 
ature  distribution  Tc0  to  that  at  the  end  of  the  noontime  heat¬ 
ing  experiment. 

can  be  modified  through  two  processes:  (1)  changes  of  the 
inelastic  collision  losses  by  the  heating  and  (2)  thermal 
expansion  along  the  geomagnetic  field. 

[h]  In  the  absence  of  D  region  absorption,  which  was  the 
experimental  situation,  the  waveficld  intensity  E()  of  the  HF 
heater  is  given  by  E0(z)  -  (3/Z)  V/m,  where  Z  =  tJzq  is  the 
altitude  z  normalized  to  a  reference  height  Zo  =  100  km. 
In  the  steady  state,  the  elevated  electron  temperature  in  a 
circularly  polarized  heater  is  given  by 

Te~Teo  +  2e2Ej/[3«(T,)mw2]  (1) 

where  u;  is  the  radian  frequency  of  the  heater,  6( Tc)  =  60  + 
£cm  is  the  electron  energy  fraction  lost  in  each  collision, 
here,  £0  =  2m/M  accounts  for  elastic  collisions  and  <$cm  = 
<5cm(Tc)  accounts  for  inelastic  collisions,  which  drain  elec¬ 
tron  energy  via  rotational  and  vibrational  excitations  of 
the  neutral  molecules  (N2  and  O2)  In  the  strong  heating 
situation,  6Cm(Te)  dominates  and  has  to  be  evaluated  self- 
consistently.  With  the  aid  of  Gurevich  [1978,  Tables  1, 2,  7, 
and  8]  and  the  electron  density  distribution  determined 
by  the  recorded  ionograms,  the  spatial  distribution  of  Tc 
is  evaluated. 

4.1.  Processes  Involved  in  Noontime  Modifications 

[15]  Presented  in  Figure  5  is  a  comparison  of  the  electron 
temperature  distribution  at  the  end  of  the  heating  experiment 
to  the  unperturbed  distribution  T^,  which  is  based  on  the 
values,  given  by  Gurevich  [1978,  Table  2].  As  shown,  the 
heater  causes  a  significant  increase  in  the  electron  temper¬ 
ature,  in  particular,  in  the  region  from  90  to  120  km  where 
the  electron  density  drops  by  as  much  as  10%. 

[16]  Electron  heating  can  also  give  rise  to  heating  of  the 
neutral  gas  in  terms  of  populating  (vibrationally)  excited  02 
and  N2.  Consequently,  the  rates  of  two  inelastic  processes 
may  be  modified  by  the  heating.  One  is  the  electron-ion 
recombination  coefficient,  which  decreases  as  the  electron 
temperature  increases;  the  second  one  is  the  eleetron-oxygen 


dissociative  attachment,  e  +  02  — ♦  O  ■+■  O  ,  in  which  the 
electron  energy  has  to  exceed  a  threshold  level  imposed 
by  the  conservation  of  energy  Normally,  this  threshold  is 
3.6  eV.  However,  this  threshold  level  decreases  as  the 
internal  energy  of  molecular  oxygen  increases  and  the  lab¬ 
oratory  experimental  results  [ Henderson  et  al ,  1969]  show 
that  this  threshold  can  decrease  to  1  eV  and  less.  Moreover, 
an  enhancement  of  the  electron  attachment  coefficient  can 
introduce  a  focusing  effect  on  the  HF  heater  [Shukla  et  al ., 
1992],  which  further  increases  the  heating  and  the  attach¬ 
ment  coefficient.  The  counterbalance  of  the  two  inelastic 
processes  on  the  electron  density  gives  a  net  effect  on  the 
electron  density  distribution. 

[17]  Therefore,  a  long  heating  period  can  change  the 
ionization  balance,  which  involves  four  processes:  (1)  photo 
ionization  by  the  solar  illumination,  (2)  artificial  ionization 
by  electrons  accelerated  by  the  strong  waveficld,  (3)  clcctron- 
heating-eaused  reduction  of  the  recombination  coefficient, 
and  (4)  enhancement  of  the  electron  attachment  coefficient 
caused  by  the  heating  of  electrons  and  oxygen  molecules. 

[18]  The  dissociative  attachment  rate  is  given  by  ia,  = 
<tNvvc,  a  product  of  the  cross  section  a,  density  Nv  of 
vibrationally  excited  molecular  oxygen,  and  the  electron 
velocity  vc,  where  vc  -  4.2  *  10  cm/s  for  a  1  eV  electron 
and  the  cross  section  o  ~~  10  21  cm2  In  the  region  from  90  to 
120  km,  1.2  x  1011  <  N(02)  <  8.2  x  1012  cm  3  [c.g., 
Gurevich ,  1978,  Table  1],  thus,  when  excited  molecular 
oxygen  constitutes  a  few  percent,  the  total  attachment  loss 
of  electrons  will  exceed  10%  of  its  density  at  the  end  of 
the  experiment.  As  the  heating  wave  propagates  upward, 
the  wavcfield  intensity  and  the  background  neutral  density 


360 


0  0  2  0  4  01  0  8  1 

(a  )  Electron  Density  (x  105  cm  3) 

360 


(b)  Electron  Temperature  (K) 

Figure  6.  Comparisons  of  (a)  electron  density  and  (b)  elec¬ 
tron  temperature  distributions  in  the  beginning  and  at  the 
end  of  the  nighttime  experiment. 
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Figure  7.  A  comparison  of  the  shifts  of  the  electron  density 
distributions  on  13  and  16  April,  caused  by  the  heating  for 
about  15  min  and  due  to  the  natural  decay  of  the  F  layer  for 
30  min,  respectively. 

decrease  with  increasing  altitude.  Thus  the  attachment  loss 
rate  decreases  with  the  altitude  When  the  positive  change  of 
the  attachment  rate  is  larger/smaller  (in  absolute  value)  than 
the  negative  change  of  the  recombination  rate,  then  the  net 
effect  can  modify  the  electron  density  distribution  with 
depiction  in  the  lower-altitude  region  and  enhancement  in 
the  higher-altitude  region,  as  shown  in  Figure  2. 

4.2.  Process  Responsible  for  the  Nighttime  Modifications 

[19]  Presented  in  Figures  6a  and  6b  are  comparisons  of 
electron  density  and  temperature  distributions  in  the  begin¬ 
ning  and  at  the  end  of  the  nighttime  heating  experiment. 


where  the  elevated  electron  temperature  distributions  in 
Figure  6b  were  again  determined  via  equation  (1).  The 
unperturbed  electron  temperature  Tc0.  from  Gurevich  [1978, 
Table  2],  is  also  plotted  in  Figure  6b  as  a  reference.  As 
shown,  the  heater  causes  considerable  electron  temperature 
elevation,  which  remains  about  the  same  through  the 
experiment.  On  the  other  hand,  the  electron  density  decreases 
continuously  in  time.  In  fact,  due  to  the  density  drop,  the 
heating  near  foF2  becomes  less  effective  in  the  later  time  of 
the  experiment.  The  results  presented  in  Figure  6  support 
the  thermal  expansion  explanation  on  the  observation.  The 
background  conditions  on  1 3  and  16  April  were  quite  similar 
as  indicated  by  the  magnetometer  measurements  and  no 
experiment  wras  performed  at  FIAARP  on  16  April  Because 
the  natural  decay  of  the  F  layer  as  the  Sun  sets  can  also  cause 
the  decay  and  upward  expansion  of  the  bottomside  electron 
density  distribution,  the  changes  of  the  distributions  after 
30  min  on  13  and  16  April  are  first  compared.  As  shown 
in  Figure  7,  the  shift  on  16  April  was  much  smaller  than 
that  on  13  April.  Next,  an  ionogram  recorded  on  16  April 
was  taken  to  compare  a  corresponding  one  (i.e.,  recorded 
at  about  the  same  time  of  the  day)  recorded  on  13  April  in 
the  experiment. 

[20]  The  true  height  profiles  and  electron  density  dis¬ 
tributions  were  also  deduced  from  the  two  ionograms  for  a 
comparison.  As  shown  in  Figures  8a  and  8b,  the  virtual 
heights  of  the  sounding  echoes  and  the  true  height  of  the 
electron  density  distribution  on  13  April  arc  much  higher 
than  the  corresponding  ones  on  16  April  The  true  height 
profiles  of  the  sounding  echoes  in  Figure  8a  were  deter¬ 
mined  by  an  algorithm  that  is  part  of  the  software  program 
SAO  Explorer  [Galkin  et  al.y  2008].  The  topside  profiles  in 
Figure  8a  arc  model  results  presented  for  reference  only,  as 
the  ionosonde  cannot  observe  the  topside  of  the  ionosphere. 
A  significant  upward  expansion  of  the  density  distribution 
on  13  April  results  in  a  considerable  density  drop  in  the 
bottomside  of  the  ionosphere.  These  two  comparisons 


Figure  8.  (a)  A  comparison  of  the  ionograms  recorded  on  13  and  16  April  and  (b)  the  deduced  electron 
density  distributions. 


6  of  7 


A03304 


KUO  ET  AL..  IONOSPHERE  UNDER-DENSE  HF  HEATING  EFFECTS 


A03304 


clearly  demonstrate  that  the  ionosphere  on  13  April  was 
lifted  mainly  by  the  heating 

5.  Discussion 

[21]  The  experimental  results  show  that  the  impacts  of  the 
under-dense  heating  by  a  powerful  HF  wave  on  the  noon¬ 
time  and  nighttime  ionosphere  are  quite  different  from  each 
other  as  well  as  different  from  that  of  the  overdense  heating 
of  the  ionosphere. 

[22]  The  overdense  heating  by  the  O  mode  HF  heater 
starts  locally  and  involves  instabilities,  such  as  parametric 
instabilities  \Fejer,  1979;  Stcnflo,  1985]  and  thermal  insta¬ 
bility  [Kuo  and  Djuth ,  1988]  and  are  frequency  dependent, 
such  as  heating  at  second  harmonic  electron  cyclotron  res¬ 
onance  [Djuth  et  al ,  2005]  Such  anomalous  heating  forms  a 
local  heat  source  which  transmits  heat  both  downward  and 
upward  through  heat  conduction  along  the  geomagnetic 
field  The  anomalous  electron  heating  causes  a  decrease  of 
the  recombination  rate,  which  changes  the  balance  between 
the  photoionization  and  the  recombination  loss  to  result  in 
an  increase  of  the  electron  density.  In  the  proeess,  a  thermal 
instability  is  also  excited  to  generate  large-scale  field- 
aligned  density  irregularities  (FAls)  [Kuo  and  Djuth ,  1988], 
which  are  responsible  for  the  observed  virtual  height  spread 
[Kuo  and  Snyder ,  2010]  The  parametrically  excited  plasma 
waves  can  also  resonantly  accelerate  electrons  to  exceed 
10  eV.  Pedersen  et  al  [2009,  2010]  suggested  those  ener¬ 
getic  electrons  to  be  the  main  cause  of  artificial  ionization 
patches  observed  in  HF  heating  experiments  at  HAARP 

[23]  On  the  other  hand,  the  under-dense  heating  covers  a 
large  altitude  region  from  the  beginning;  EIRP  of  the  heater 
is  probably  the  most  important  ingredient  Both  bottomside 
and  topside  of  the  ionosphere  can  be  perturbed  simulta¬ 
neously,  It  causes  a  larger-seale  ionospheric  modification 
which  also  lasts  longer.  The  experimental  results  show'  that 
the  heating  effect  accumulates  in  time  even  in  the  inter¬ 
mittent  heating  situation  with  a  2  min  heater-off  period 
introduced  between  two  2  min  heater-on  periods. 

[24]  In  the  noontime,  the  ionosphere  is  at  a  low'  altitude, 
thus,  collisions  are  dominated  by  the  inelastic  processes.  On 
the  other  hand,  elastic  collisions  dominate  in  the  nighttime. 
Therefore,  the  electron  temperature  elevation  in  noontime 
under-dense  heating  decreases  with  an  increase  in  altitude 
and  the  trend  is  reversed  in  the  nighttime  The  electron 
temperature  increases  significantly  in  both  experimental  time 
periods.  The  results  of  the  noontime  experiment  show'  that  the 
electron  density  deereases/inereases  in  time  in  the  region 
below/above  a  critical  height  at  about  140  km,  manifesting  a 
change  of  the  balance  between  photoionization  and  recom¬ 
bination  and  attachment  losses  by  the  heating.  The  results  of 
the  nighttime  experiment  show'  that  the  ionosphere  was  lifted 
by  30  to  50  km  through  continuously  upward  expansion.  In 
the  future  work,  a  numerical  simulation  effort,  similar  to  that 
[Eliasson  and  Stenflo ,  2010]  for  the  phenomena  explored  in 
the  overdense  heating  case,  is  also  called  for 
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